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Parity-time-symmetric (PT-symmetric) optical waveguide couplers have become a key component 
for integrated optics. They offer new possibilities for fast, ultracompact, configurable, all-optical sig¬ 
nal processing. Here, we study nonlinear properties of finite-size multimode PT-symmetric couplers 
and report a peculiar type of dispersion relation in couplers made of more than one dimer. More¬ 
over, we predict a VT transition triggered by nonlinearity in these structures, and we demonstrate 
that with the increase of the number of dimers in the system, the transition threshold decreases and 
converges to the value corresponding to an infinite array. Finally, we present a variety of periodic 
intensity patterns that can be formed in these couplers depending on the initial excitation. 

PACS numbers: 42.65.Wi, 42.79.Gn, 42.65.Sf 


Quantum mechanical observables correspond to oper¬ 
ators with a real spectrum of eigenvalues [T]. Hermitian 
operators are the most well-known class of such opera¬ 
tors. However, a more general class of non-Hermitian 
operators that possess a real spectrum was described by 
Bender et al. 013 ] in 1998. These parity-time-symmetric 
(T^T-symmetric) operators commute with a subsequent 
action of the parity operator V and the time inversion 
operator T. Consequently, T^T-symmetric Hamiltonians 
have real spectra and describe physical phenomena, pro¬ 
vided that the complex potentials V fulfill the condition 
y(r) = ]/*(—r). Optical systems, where the light prop¬ 
agation is described by the nonlinear Schrodinger equa¬ 
tion, offer an efficient platform for implementation of VT 
symmetry. The complex dielectric permittivity distribu¬ 
tion e(r) plays the role of potential and the imaginary 
part of the permittivity corresponds to gain or loss ISIS]. 

The optical VT symmetry in the linear regime was 
experimentally demonstrated in waveguide couplers [6], 
which have a critical importance for the future ultracom¬ 
pact, fast, and configurable all-optical switching. In the 
nonlinear regime, optical couplers with gain and loss were 
studied [7] , where suppression of time reversal [8] and uni¬ 
directionality [9] were demonstrated. Phenomena related 
to nonlinear VT symmetry were also investigated in pe¬ 
riodic systems, where solitons HqHISI, breathers [16] and 
their stability iniiTi] were analyzed. 

For fixed geometrical parameters, T^T-symmetric sys¬ 
tems can be in the full or the broken PT-symmetric 
regime depending on the ratio between the imaginary 
eiM and the real part eRE of the permittivity modulation 
depth (eiivi/eRE)- For low values of this ratio, the sys¬ 
tem is in the full T^T-symmetric regime and has purely 
real eigenvalues mm- By increasing the eiu /^re ratio, 
the VT transition threshold is reached, and the system 
transforms into the broken T^T-symmetric regime, where 
a pair of modes (one with gain and one with loss) has 
complex conjugate effective indices. In the linear regime, 
the €im/^re ratio is changed by varying ejM, responsi¬ 


ble for gain and loss in the system. However, changes of 
the real part of the permittivity also influence the ratio 
of eiivi/eRE- The amplitude of eRE in nonlinear systems 
can be controlled by varying the incident light intensity. 
Recently, a nonlinearity-triggered transition from the full 
to the broken T^T-symmetric regime (VT transition) was 
reported in an infinite periodic array of waveguides de¬ 
scribed by cosine-like permittivity distribution [20] . 

In this paper, we study the nonlinear, T^T-symmetric, 
ultracompact couplers consisting of a few cosine-like 
waveguide dimers that can be readily integrated on a 
chip. Figure shows the three-dimer geometry under 
investigation. We report a new type of linear disper¬ 
sion relation for a multimode PT-symmetric coupler [21] 
built of more than one dimer. Additionally, we predict a 
VT transition in such systems and study the dependence 
of the transition threshold on the number of dimers in the 
coupler. Peculiar mode interference patterns are reported 
for the coupler built of three dimers. 

Light propagation in one-dimensional (ID) VT- 
symmetric structures with cubic nonlinearity is studied 
using a scalar wave equation for the electric field E(x, z): 

+ k^€{x) + k^a{x)\Ef] E = 0, (1) 

where ko = 27r/A is the free-space wavevector, A is the 



FIG. 1. Geometry of the studied PT-symmetric coupler with 
its parameters. The size (period) of the cosine-like dimers is 
denoted by F. cb denotes the background relative permittiv¬ 
ity; cre and eiM denote the modulation amplitude of the real 
and the imaginary part of relative permittivity, respectively. 
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free-space wavelength of light, and the operator = 
d‘^ jdx^ + d‘^ jdz^ denotes the 2D Laplacian. The light 
propagates along the 2 ;-direction, and both the structure 
and the field distributions are assumed to be invariant 
along the ^-direction. The linear complex relative per¬ 
mittivity distribution is described by e{x) = + Ae(x), 

where es denotes the background relative permittiv¬ 
ity and Ae(x) describes the linear complex permittivity 
modulation depth. The Kerr-type nonlinearity strength 
is quantified by the parameter a. The nonlinearity a{x) 
is nonzero only inside of the waveguides . 

Using the slowly varying envelope approximation 
E{x,z) = '0(x, Eq. § is transformed into 

the (1 -h 1)D nonlinear Schrodinger equation 


dilj 

dz 




^^+^0 (Ae + a|^H 




( 2 ) 


In order to analyze the nonlinear dynamics of the 
light propagation in T^T-symmetric waveguides, we solve 
Eq. § using the split-step Eourier method [22 [23] . 

Previously, nonlinear VT transition was theoretically 
studied in an infinite periodic cosine-like waveguide ar¬ 
ray [20]. Here, we analyze simpler, finite-size systems 
of waveguides, where we find a rich variety of nonlinear 
phenomena. Let us consider a finite-size array built of 
three dimers (shown in Eig. described by 

Ae{x) = eRE cos^(Ux) + ieiu sin(2Ux), (3) 


where U = tt/T, T denotes the full width of the dimer, 
and the structure parameters are the following: T = 
3 /rm, eRE = 0.05, and a = 10“^^ m^/V^ at A = 0.63 /rm. 

Let’s first consider a linear system [a{x) = 0]. The 
field profiles that propagate with an effective index Ueff 
without changing their shape are sought in the form 
E{x,z) = This allows ns to transform 

Eq. § into an eigenvalue problem for the modes of 
the waveguide described by e{x). The eigenvalue prob¬ 
lem is solved using a finite-difference method resulting 
in the dispersion relations [shown in Eigs. [^a), (b)] and 
the corresponding field profiles of the eigenmodes [shown 
in Eigs. I^c), (d)]. We use these dispersion diagrams to 
locate the VT transition threshold. 

Eor the parameters of the dimers chosen here, the 
waveguide supports more than one mode. Eor low values 
of eiM, there exist four modes with purely real effective 
indices. With the increase of ejM, two new modes emerge 
at eiM ~ 0.007 and 0.34. Meanwhile, the pairs of modes 
with real effective indices undergo a transition to the bro¬ 
ken T^T-symmetric regime. Eor high values of ejM, the 
higher-order conjugated mode pairs vanish, leaving only 
the first mode pair from which one mode experiences gain 
and the other experiences loss during the propagation. 

In a single multimode dimer [Eig. |Sl[a) in supple¬ 
mentary material], the dispersion curves of higher-order 
mode pairs appear below the dispersion curves of the two 
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FIG. 2. (a), (b) Dispersion diagrams for a T^T-symmetric cou¬ 
pler composed of three dimers: (a) the real and (b) the imag¬ 
inary part of the effective index of the modes as a function of 
the imaginary part of permittivity eiM in the linear case. The 
black dashed line denotes the refractive index of a background 
medium yTT- (c), (d) Absolute value of the field distribu¬ 
tions |0(x)| of (c) the lossless modes [mode 1 (fundamental)— 
red, mode 2—green, mode 5 (lowest effective index)—blue], 
(d) the coupled modes with loss (green) and gain (red) at 
eiM = 0.027. Black curves indicate 3^e{Ae} (solid) and 
^m{Ae} (dashed). 


first modes (see, Eig. S2 in supplementary material and 
Eigs. 2, 5, and 7 in Ref. [21]). On the contrary, in the 
case of an array built of multiple dimers, the dispersion 
curves of the mode pairs with purely real effective indices 
do not lay above each other but inside (dispersion curves 
of the lower-order modes are enclosed by these of the 
higher-order mode pairs), as it is seen in Eig. Above 
the VT transition, the curve corresponding to 5Re{neff} 
of the two lower-order modes intersects with the disper¬ 
sion curve of one of the lossless higher-order modes [see 
e.g., the point labeled A in Eig. [^a) corresponding to 
eiM = 0.0314]. At these points, three modes with the 
same value of 5Re{neff} exist simultaneously in the dimer 
array. One of them experiences gain, one loss, and one 
propagates without change of its amplitude. 

The dispersion relations presented in Eig. [^provide the 
value of eiM where the first (for lowest value of cim) VT 
threshold occurs in the system. Eor the array of three 
dimers, this value is cjm = 0.0268. Below the thresh¬ 
old we find five symmetric (with respect to the center 
of the array x = 0) lossles modes with real effective in¬ 
dices. Above the threshold, in addition to the three loss¬ 
less modes, we find two modes: one for which the field 
profile is shifted slightly to the gain region (gain mode) 
and one shifted to the loss region (loss mode). The cor¬ 
responding field profiles are presented in Eigs. He), (d). 
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FIG. 3. Nonlinear dynamics of light propagating in an array of 
three dimers described by Eq. § and shown in Fig. 0 (a), (c) 
The intensity distribution I{x,z), (b), (d) the evolution of 
the center of mass. eiM is equal to 0.027. The input is the 
linear gain mode with the power density Pq equal to: (a), (b) 
5 • 10® W/m and (c), (d) 10® W/m. 


In addition to the coupler built of three dimers studied 
here, we have analyzed the VT transition in structures 
made of different numbers of dimers. Here, we summarize 
the values of the VT transition threshold obtained: 
for a single dimer = 0.0305 [24]; for two dimers: 
0.0278 (see supplementary materials for more details), 
and for three dimers: 0.0268. Additional calculations 
show that for the case of four dimers, this value is equal 
to 0.0262. This suggests two conclusions. Firstly, the 
threshold value for a finite-size array of dimers is higher 
than the threshold for the infinite array of dimers de¬ 
scribed by Eq. which was reported in Refs. nniiio] 
to be = O.Scre = 0.025. This can be understood 
considering a simple model based on the fact that light 
concentrates in the regions with high permittivity. Then, 
the rate with which the overlap integral between a test 
field profile and the permittivity profile decreases, while 
shifting the test field profile along the x-direction quan¬ 
tifies how difficult it is to break the symmetry in the 
system and induce the VT transition. The higher the 
number of dimers, the slower the change of this overlap 
integral is and, consequently, the easier it is to induce 
the T^T-symmetry-breaking transition. The second con¬ 
clusion states that with the increase of the number of 
dimers, the threshold value converges to the value cor¬ 
responding to the infinite array as expected because the 
structure becomes more similar to the infinite array. 

Figure presents the evolution of the light propagat¬ 


ing in the nonlinear coupler built of three T^T-symmetric 
dimers described by Eq. (§. Initially, the system is in 
the broken T^T-symmetric regime, and the light is in¬ 
jected into the gain mode of the system. The inten¬ 
sity of this mode grows rapidly, causing the nonlinear 
increase of the real part of permittivity, which trans¬ 
forms the system to the full T^T-symmetric regime. In 
this regime, the light is attracted toward the center of 
each dimer (where the permittivity is the highest), and 
due to the transverse momentum, light crosses the cen¬ 
ter and is mostly located in the loss region m- Conse¬ 
quently, the light intensity and the nonlinear modifica¬ 
tion of permittivity decrease bringing the system back to 
the broken T^T-symmetric regime. This cycle repeats as 
shown in Fig. [^a) and the evolution of the x-coordinate 
of the center of mass (COM) of the intensity distri¬ 
bution COM{z) = xl{x, z)dx/ I{x, z)dx pre¬ 
sented in Fig.j^b). For the input beam energy chosen in 
Figs.j^a) and (b), the total power in the system increases 
25 times during one cycle, and the period of the cycle is 
approximately equal to 1 mm. Lowering the input power 
results in an increase of the propagation length required 
for the power to grow and initiate the oscillations, but 
the oscillation period and the maximum power remains 
approximately the same. On the contrary, an increase of 
the input power results in a decrease of the oscillation 
period (down to 0.5 mm) and an increase (1.5 times) of 
the peak power [see Figs. ic), (d)]. 

Figure shows nonlinear propagation of light in the 
multimode system studied in Figs.[^and[^ but this time 
instead of exciting a single (gain) mode, we excite a mix¬ 
ture of modes. Figure [^a) shows the result of excitation 
with a wide Gaussian beam illuminating all three dimers. 
This beam excites all the modes of the waveguide (except 
the antisymmetric mode 2), and these modes interfere 
during the propagation. This interference results in a 
doubly periodic light distribution pattern. The high fre¬ 
quency oscillations result from the interference of all the 
modes that have different effective indices. Additionally, 
we observe low frequency oscillations, that are the result 
of the nonlinear VT transition that couples light between 
the gain and the loss modes. The resulting pattern resem¬ 
bles the triangular lattice. A similar interference pattern 
is observed when a mixture of the fundamental mode and 
the gain mode are excited [see Fig. [^b)]. 

On the contrary, the excitation of the mixture of 
mode 2 and the gain mode results in a different in¬ 
terference pattern, as shown in Fig. [^c). Mode 2 is 
mainly localized in the side waveguides [see green curves 
in Fig. I^a)]; therefore, the high frequency interference 
pattern is observed in these two waveguides only. The 
low frequency pattern attributed to the energy exchange 
between the gain and the loss mode is spread over all 
three waveguide dimers similar to earlier cases. Here, 
the pattern is less regular than before and consists of an 
elongated maximum in the middle waveguide and a few 
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FIG. 4. The intensity distributions I{x^z) showing the non¬ 
linear dynamics of light propagating in an array built of three 
dimers with the same parameters as these presented in Fig. 
The input field is in the form of (a) a wide Gaussian beam 
illuminating three dimers, (b) the sum of mode 1 and gain 
mode, (c) the sum of mode 2 and gain mode (for mode num¬ 
bers and profiles see Fig. 0. In all the cases, the excitation 
power density is Po = 5 • 10^ W/m. 


shorter maxima in the side waveguides. 

In conclusion, we have studied nonlinear dynamics 
of light propagation in finite-size parity-time-symmetric 
couplers built of three dimers, and reported a new type 
of linear dispersion diagram neff(eiM) for multimode sys¬ 
tems built of several dimers, where the dispersion curves 
of higher-order mode pairs are enclosed in the disper¬ 
sion curves of lower-order mode pairs. This contrasts 
with the case of multimode single dimers, where the dis¬ 
persion curves of higher-order modes appear below the 
curves corresponding to lower-order modes. Moreover, 
the nonlinearly triggered transition from the full to the 
broken PT-symmetric regime can be observed in such a 
system. The threshold for this VT transition was found 
to be higher in finite-size systems than for infinite peri¬ 
odic potential. Finally, we predicted a variety of periodic 
patterns can be formed in the multimode coupler by con¬ 
trolling the initial excitation conditions. 
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Supplementary materials 


In the main text of the manuscript, we have studied the 
nonlinear T^T-symmetric coupler built of three dimers. 
Here, we compare its dispersion relations and the nonlin¬ 
ear dynamics of light propagation with simpler structures 
built of one or two dimers described by Eq. (§ in the 
main text. The geometries of the structures described in 
these supplementary materials are presented in Fig. |S1[ 



FIG. SI. Geometry of the PT-symmetric couplers built of 
(a) one and (b) two cosine-like dimers. The period of the 
cosine-like dimers is denoted by F. cb denotes the background 
relative permittivity; cre and eiM denote the modulation am¬ 
plitude of the real and the imaginary part of relative permit¬ 
tivity, respectively. 


Linear dispersion diagrams for a single multimode dimer 


dimer are F = 3 /im, = 2, and cre = 0.3. The modu¬ 
lation depth of the real part of permittivity is chosen in 
such a way that the dimer supports more than one pair of 
modes. In this case, the dispersion curves of higher-order 
modes appear above the dispersion curves of the lowest- 
order modes, as it was shown in Ref. m- In contrast 
to the dispersion relation of a single multimode dimer 
presented in Fig. the dispersion relations of a mul¬ 
timode coupler that is built of more than one dimer, as 
shown in Fig. |^in the main text and in Fig. |S3[ present 
a different type of behavior. In this system, the disper¬ 
sion curves of the higher-order modes appear inside the 
curves corresponding to the lower-order modes. 


Coupler built of two dimers 


In the main text, we have presented the description of 
a nonlinear coupler described by a cosine-like permittiv¬ 
ity profile [Eq. Q in the main text] built of three VT- 
symmetric dimers. Here, we complete the description by 
studying an array composed of two dimers, whose ge¬ 
ometry is presented in Fig. [Hb). The linear disper¬ 
sion relations of the system of two dimers, presented 


in Fig. S3 [a), (b), show qualitatively similar behavior to 
the coupler built of three dimers. Here, the novel type 


In the main text, we have reported a new type of 
dispersion diagrams for the multimode couplers built of 
more than one dimer. Here, we compare these new dis¬ 
persion diagrams [shown in Fig. [^a), (b) in the main 
text] with the dispersion relations of a multimode waveg¬ 
uide composed of a single T^T-symmetric dimer. Fig¬ 
ure presents the dispersion relations neff(eiM) for a 
single dimer shown in Fig. |Sl[a). The parameters of the 



FIG. S2. Linear dispersion relations for a single cosine-like 
dimer described by Eq. § and shown in Fig. |Sl[ a). The (a) 


real part and (b) imaginary part of the effective index are 
shown as a function of the imaginary part of the permittivity 
modulation depth eiM- Dimer parameters (given in the text) 
are chosen in such a way that the dimer supports multiple 
modes. The black dashed line denotes the refractive index of 
a background medium yTT. 
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FIG. S3, (a), (b) Dispersion relations for a T^T-symmetric 
coupler composed of two dimers; (a) the real and (b) the imag¬ 
inary part of the effective index of the modes as a function of 
the imaginary part of permittivity cim in the linear case. The 
black dashed line denotes the refractive index of a background 
medium yTT. (c), (d) Absolute value of the field distribu¬ 
tions \<f>{x) \ of (c) the lossless modes [mode 1 (fundamental)— 
red, mode 4 (lowest effective index)—black], (d) the coupled 
modes with loss (blue) and gain (green) at eiM = 0.028. Black 
curves indicate 5Re{Ae} (solid) and ^m{Ae} (dashed). 
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FIG. S4. The intensity distributions I{x, z) showing the non¬ 
linear dynamics of light propagating in an array built of two 
dimers described by Eq. ©• The geometry of the structure is 
shown in Fig. SI and in the inset of supblot (a). The parame¬ 
ters are: T = 3 /im, cre = 0.05, a = 10“^® m^/V^ and (a), (b) 
eiM = 0.028, (c), (d) eim = 0.030. Excitation power density 
Po is equal to: (a), (c) 10^ W/m and (b), (d) 10® W/m. The 
propagation distance shown in (a) is three times longer than 
in (b)-(d). 


permittivity modulation depth increases to 1.8 • 10“^. In 
the case of two dimers, the effect of the increase of the 
amplitude of the imaginary part of permittivity is also 
studied. Figures [S^c), (d) present the light propaga¬ 
tion for the same two input powers as in Figs. |S4[a), (b) 
but in the system located higher above the VT tran¬ 
sition. In this case, the oscillation period decreased to 
0.25 mm. Contrary to the case of lower amplitude of the 
imaginary part, the increase of the initial power for the 
case presented in Figs. |S4[c) and (d) does not modify 
the oscillation period, and it is the same for both power 
levels. Only the initial propagation distance required to 
initiate the oscillations (from the input plane 2 : = 0 to 
the first intensity maximum) is reduced with the increase 
of the excitation power because the same intensity level 
is reached earlier by the system, where more power was 
injected. 

Finally, we have studied the behavior of the coupler 
built of two T^T-symmetric dimers when more than one 
mode of the system is excited. The field profiles of 
the linear modes supported by the system are presented 
in Fig.|S3[c), (d). Figure [^presents the evolution of the 
field in the case when the gain mode and the fundamental 
mode are excited simultaneously. The interference pat¬ 
tern obtained here resembles these for the coupler built 
of three waveguides presented in Fig. in the main text. 
However, for the simpler system composed of two waveg¬ 
uides, the interference pattern is less complex than for 
the coupler built of three T^T-symmetric waveguides. 


of dispersion diagrams is also observed, where the dis¬ 
persion curves of the higher-order mode pairs enclose the 
curves corresponding to lower-order mode pairs. This 
confirms our conclusions that the dispersion curves of a 
multimode T^T-symmetric coupler built of a single dimer 
are qualitatively different from the curves of couplers 
built of more than one dimer. 

Similar to the case of a coupler built of three dimers, 
here we also observe the nonlinearly-triggered-T^T tran¬ 
sition from the full to the broken T^T-symmetric regime. 
Linear modal studies of this system show that this tran¬ 
sition occurs at ejM = 0.0278. The nonlinear propagation 
of light in the system just above the VT transition at low 
and high input power levels is shown in Figs. [S4|^a) and 
(b), respectively. At the low power level, the oscillation 
period is approximately equal to 1.4 mm, and the peak 
power corresponds to the nonlinear permittivity modu¬ 
lation depth of the order of 6 • 10“^. For high power, the 
period decreases to 0.4 mm, and the maximum nonlinear 
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FIG. S5. The intensity distributions I{x,z) showing the non¬ 
linear dynamics of light propagating in an array built of 
two dimers with the same parameters as these presented in 
Fig. [S4|b). The input held is in the form of the sum of 
mode 1 and gain mode. The ex citat ion power density is 
Po = 10® W/m [the same as in Fig. S4[b)]. 




















